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We propose a new quantum algorithm for finding eigenvalues of non-Hermitian matrices. We show how to 
construct, through interactions in a quantum system and partial measurements, a non-Hermitian matrix and ob- 
tain its eigenvalues and eigenvectors. This proposal combines ideas of frequent partial measurement, measured 
quantum Fourier transform, and quantum state tomography. Similar to the conventional phase estimation algo- 
rithm, our algorithm provides a speed-up over its classical counterpart. Moreover, in our approach the controlled 
unitary operation is implemented only once. In addition, the algorithm provides a high success probability. 

PACS numbers: 03.67.Ac, 03.67.Lx 
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Introduction - One of the most important tasks for a 
quantum computer is to efficiently obtain eigenvalues and 
eigenvectors of high-dimensional matrices. It has been 
suggested [JJ that the quantum Phase Estimation Algo- 
rithm (peaTiI can be used to obtain eigenvalues of a Her- 
mitian matrix or Hamiltonian. For a quantum system with a 
Hamiltonian H, a phase factor, which encodes the informa- 
tion of eigenvalues of H, is generated via the unitary evolu- 



tion U 
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By evaluating this phase, we can obtain the 



eigenvalues of H. The standard PEA consists of four steps: 
preparing an initial approximate eigenstate of the Hamiltonian 
H, implementing unitary evolution operation, performing the 
inverse quantum Fourier transformation (QFT), and measur- 
ing binary digits of the index qubits. 

The PEA is at the heart of a variety of quantum algorithms, 
including Shor's factoring algorithm [3]. A number of ap- 
plications of PEA have been developed, including generating 
eigenstates associated with an operator |4], evaluating eigen- 
values of differential operators [5^, and it has been generalized 
using adaptive measurement theory to achieve a quantum- 
enhanced measurement precision at the Heisenberg limit [6]. 
The PEA has been applied in quantum chemistry to obtain 
eigenenergies of molecular systems 10, Ht]. This application 
has been demonstrated in a recent experiment (9] . 

PEA is one of the few quantum algorithms that provide 
speed-up over their classical counterparts. It is only designed 
for finding eigenvalues of a Hermitian matrix. However, in 
modern engineering and science, people are often more inter- 
ested in eigen-problems of non-Hermitian matrices. In this pa- 
per, we propose a measurement-based phase estimation algo- 
rithm (MPEA) to evaluate eigenvalues of non-Hermitian ma- 
trices. Our proposal concatenates ideas in conventional PEA, 
partial measurement, and current techniques in one-qubit state 
tomography. The method in this proposal can be used to de- 
sign useful quantum algorithms apart from those based on the 
standard unitary circuit model. Like the conventional PEA, 
we expect our MPEA to play a key role in discovering new 



quantum algorithms. 

Constructing non-Hermitian matrices. -Now we describe 
how to construct non-Hermitian matrices for a quantum sys- 
tem. A bipartite system, composed of subsystems A and B, 
evolves under H = Ha + Hb + Hab, where Ha (b) is 
the Hamiltonian of subsystem A{E) and Hab is their inter- 
action. If subsystem A is subject to a projective measurement 
M — \ip^){ip^ \ applied during a time interval r, where \ip^) 
is a pure state of subsystem A, this is equivalent to driving sub- 
system B with an evolution matrix Vb{t) = {ip^\e^^^'^\ip^). 
This evolution matrix is in general neither unitary nor Her- 
mitian. We can construct such a matrix Vb by choosing an 
appropriate basis for the partial measurements. 

The Hamiltonian H of the whole quantum system can be 
spanned as: H = w'\\h eigenenergies Ej. 

The eigenvector \Ej) can be spanned in terms of tensor prod- 
ucts of basis vectors {IV'fe )} and of Hilbert spaces of 
subsystems A and B, which are of dimensions Da and Db, 
respectively, and D = Da 'x Db- Using basis for A and B, 
one can expand \E,) - Y.k=i Y.r=i fl,M) ® l^f )' 
the evolution matrix on subsystem B, after a measurement on 
A during a time interval r, becomes 
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((/? iV'fe )• We can thus construct any non-Hermitian matrix 
that can be written in this form on a quantum system for any 
given arbitrary /^^ and |v5"*). 

More generally, we can construct a different evolution 
matrix by performing measurements on subsystem A with 
different time intervals and/or different measurement bases. 
For example, by sequentially performing partial measure- 
ments with time intervals t\, T2, ts, an evolution matrix 
VB(ri,T2,T3) = Vb{tj,)Vb{t2)Vb{t\) is constructed. We 
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can also combine unitary evolution matrices with the non- 
Hermitian transformations on subsystem B to construct some 
desired evolution matrices. 

Measurement-based quantum phase estimation algorithm.- 
We now present a MPEA for evaluating the eigenvalues and 
eigenvectors of a non-Hermitian matrix Vb(t). There are 
three primary steps in the MPEA: (1) preparing an initial 
state; (2) constructing controlled successive binary powers of 
the evolution matrix Vb{t); and (3) retrieving the eigenvalue 
of the evolution matrix Vb(t) using either quantum state to- 
mog raphy (QST) | or a measured quantum Fourier trans- 
form (mQFT) combined with projective measurements. 

We now prepare the system in a separable state = 
Wa){'Pa\ ® PB' with A in a pure state |(^^) and B in an 
arbitrary (mixed) state pg. Constructing the controlled evo- 
lution matrix of Vb(t) on B is achieved by implementing 
the Controlled Unitary transformation (C-U), for the whole 
quantum system and performing the periodic measurement 
M — |'P^)(</5^| on A with period r. In order to resolve the 
phase factor up to n binary digits using the inverse QFT, one 
has to construct a series of evolution matrices of successive bi- 
nary powers, from (n — 1) to 0. In the MPEA, this is done by 
implementing the C-U operation for the whole system once 
and performing a series of 2'^ periodic measurements in time 
intervals r for a number of successive binary powers of k. 
Then we can obtain a series of transformation matrices in bi- 
nary powers of Vb(t), [yB(T)]^ , A; = (n— 1), (n— 2), • • • ,0. 

It can be shown that after m successful partial measure- 
ments on A, the evolution operation on the Hilbert space of 
subsystem B is {Vb (t))'" and the state of B evolves to ifTTIl 



'PB[vUr)r/Pn 
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where ^Ttb{[Vb {t)]"'pb [V^ (t)]™} is the probabiUty of 
finding subsystem A still in its initial state \(pj^) after m mea- 
surements. The evolution matrix Vb(t) can be spanned as 
^s(t) = J2k ^k\uk){vk\, where \uk) and {vk\ are the right- 
and left-eigenvectors and Xk is the corresponding eigenvalue 
satisfying < |Aa;| < 1 for arbitrary k Ollll . In the large 
TO limit, the operator [Vb (t)]™ is dominated by a single term 
^maxl^max)(winax|, providcd the largest eigenvalue A^ax is 
unique, discrete, and nondegenerate. In the limit of large m 
and finite r, p^' (to) tends to a pure state, independent of the 
initial (mixed) state of system B. The final state of p^'' (to) is 



dominated by 
ability — 



), and this outcome is found with prob- 



<) {v„ 



The 



State of subsystem B evolves to |umax), after performing a 
number of operations of Vb{t). Then we can evaluate A,„ax 
by resolving the phase of the state. 

The circuit for the MPEA is shown in Fig. 1. In Fig. 1(a), 
we show the circuit for fcth partial measurement with period 
T, W{k), where a measurement M = \iPa){'Pa\ is performed 
2'' times in the period r, as the whole system evolves under 
the controlled unitary operation. Fig. 1(6) shows the MPEA 
circuit; from top to bottom: three quantum registers are pre- 
pared: an index register, a target register and an interacting 



register The index register is a single qubit, which is used as 
control qubit and to obtain the final results; the target register 
is used to represent the state of subsystem B; the interacting 
register is used to represent the state of subsystem A, and to 
perform measurements. 

The initial state of the MPEA circuit is prepared in state 
(|0) + |1))((0| + (1|) (E)Pb'» \'Pa){'PaI after performing the 
(n — l)th periodic measurements W{n — 1) on the interact- 
ing register, the state of the system is transformed to (|0) + 

\l)[VB{r)r'')m + {l\K{r)r-') <^_Pb ® \Va){VaI 
The dominating term is (|0) + (Amax)^" |l))|umax)|<PA>- 
There are two approaches to resolve Amax at this stage. In 
QST, we prepare a large number of identical copies of the 
state on the index qubit, we can obtain A,nax by determin- 
ing the index qubit state. In the second approach, we com- 
bine the mQFT technique and projective measurement. The 
phases which encode the eigenvalues of Vb{t) in general are 
complex numbers, the QFT can only help resolving the real 
part of the phase. In this approach, we perform measure- 
ments on subsystem A for n series in sequence of W{n — 
1), W{n — 2), • • • , M^(0), then we obtain single qubit states 

(|0) + (A,„ax)^'""'' |1)), (|0) + (A„,ax)2'""'' |1)), • • • , (|0) + 

Amax 1 1)) on the index qubit. The real part of the phase factor 
of Amax can be obtained by performing mQFT. This proce- 
dure is shown in Fig. 1, where the circuit in the dashed square 
in Fig. 1(6) is replaced by circuits in the dashed squares shown 
in Fig. 1(c) sequentially to obtain n binary digits of the phase. 
The technique for resolving the imaginary part of the phase 
factor will be discussed in the following sections. Next we 
give details of obtaining all eigenvalues of Vb(t). 

Obtaining complex eigenvalues through state tomography.- 
QST can fuUy characterize the quantum state of a particle 
or particles through a series of measurements in different 
bases |13]. In MPEA, after performing a large number, to, 
of partial measurements, the state of the index qubit is dom- 
inated by IV-ind) = |0) + (Amax)™ 1 1)- In general Amax can 
be written as Amax — e**^ — pi{a+ib') ^ obtain Amax 

by retrieving the phase factor ip = {a + ib). Using QST, 
we prepare a large number of identical copies of the sin- 
gle qubit state IV^m) = |0) + e""'''e*''""|l). To determine 
this state, we can perform a projective measurement in the 
basis |1)(1| to obtain the probability of finding the state as 
|1), thus obtaining the value of 6. Performing a 7r/2 rotation 
around the x-axis and a measurement in the basis of the Pauli 
matrix on the index qubit, we can obtain the observable 



and thus obtain the value of a. 



The measurement errors of QST, from counting statistics, 
obey the central limit theorem. To obtain more accurate re- 
sults, we have to prepare a larger ensemble of the single qubit 
states. 

Obtaining complex eigenvalues through mQFT and projec- 
tive measurements.- In the MPEA circuit shown in Fig. 1, 
we resolve the real part of the phase factor using mQFT. 
The mQFT technique implements QFT using only a single 
qubit [10.1 . It uses the fact that gates within the Fourier trans- 
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FIG. 1: Quantum circuit for measurement-based phase estimation 
algorithm (MPEA). Part (a) shows the circuit for the "fcth" partial 
measurement Wik) with period t; part (6) shows the circuit for the 
MPEA. From top to bottom, an index register, a target register and 
an interacting register are prepared in the states l/\/2(|0) + |1)), 
Pg, and \'~PaI, respectively. The index register is a single qubit used 
as control qubit and to perform mQFT; the target register is used 
to represent the state of subsystem S; and the interacting register is 
used to represent the state of subsystem A, which interacts with B. 
The circuit in the dashed square is used for obtaining the nth binary 
digit of the real part of the phase factor In part (c), a sequence of 
circuits inside the dashed squares is used to replace the circuit in the 
dashed square in part (6), in order to resolve from the (n — l)th to the 
1st binary digits of the real part of the phase factor Here, H is the 
Hadamard gate, and _R„ is a single-qubit rotation. 



form are applied sequentially on qubits. This modification of 
the QFT algorithm preserves the probabilities of all measure- 
ments 1 12]. We employ a projective measurement to resolve 
the imaginary part of the phase factor. 

In the MPEA, each time when we perform the peri- 
odic partial measurements Wik), we obtain a state |0) + 
(Amax)^'^ll) on the index qubit. We rewrite this state as 
IV'ind) = ^(|0) +e^('^+^^)|l)) = -i=(|0) + e-Vll)). 

Let (1/V2)2 + (e-VV2)2 = r^, then we have IV'ind) = 
— ( 1 0) +e^''e*'' 1 1) ) . We now prepare a large number of iden- 
tical copies of iV'ind) by running the MPEA for many times. 
Then, performing a projective measurement |0) (0| on the en- 
semble of qubits, we can evaluate the probability for observ- 
ing the qubit in state |0), and therefore obtain r and h. 

After we obtain the value of 6, we run MPEA again and 
perform a single qubit operation on the index qubit such that 



Then we apply the mQFT technique to resolve the real part a 
of the phase factor a + ih. We therefore obtain the eigenvalue 
gj(a+jb) jjjg non-Hermitian matrix Vb{t)- 

In the above procedure, the state on the target register is 
driven to the eigenstate with the largest eigenvalue of Vb (t) 
by performing periodic partial measurements. The target state 
can evolve to the eigenstate with the second largest eigen- 
value of Vb(j) by projecting out |wniax) in preparing the ini- 
tial state. If we prepare the initial state on the target state as 



I Pb I Umax)), and feed it into the MPEA 



circuit, we obtain the second largest eigenvalue of yB(r). All 
eigenvalues of Vb (t) can be obtained by repeating this proce- 
dure. 

The MPEA resolves the nth binary digit of the phase factor 
first, and the partial measurement on A is performed in se- 
quence of 2"^^, 2"^^ to "2? times. This procedure provides 
high fidelity for the state of subsystem B since each measure- 
ment drives the state of B closer to the eigenstate of Vb(j\ 
The success probability of the MPEA is the product of the 
fidelity of the state and the probability of successful measure- 
ments on subsystem A. Since this fidelity is close to one, the 
success probability only depends on the probability of suc- 
cessful measurements on A. 

An example - Let us now use a simple model to show 
how MPEA works. We consider a quantum system con- 
sisting of two subsystems A and B, where B contains two 
non-interacting spin qubits, B\ and B2; and A is a photon. 
The whole system can be described by the Jaynes-Cummings 
Hamiltonian [15] 

(3) 

where b (6^) is a bosonic annihilation (creation) operator of 
the photons. Consider the case wq = wi, and perform partial 
measurements in the basis of a single photon state | (/3^) = 1 1) . 
We thus have VbIt) = diag{l, [3 + 2 cos(V10t J)]/5 x 
e-2™or ^ ^Qg^ygr J)e-™o^, cos(v^t J)}, in the ordered ba- 
sis {|1, s), |1, i+), |1, io), |1, L)}, where |s) = 73(101)- 

|10)), \t+) = 111), |<o) = ^(|01) + |10)) and \t_) = |00). 
Let r = 1/2 and wo ^ wi = J = 1, then, the evolution ma- 
trix (r) becomes (r) = diag{l, [3 + 2cos(\/T0/2)]/5x 
e~\ cos(\/6/2)e~*/2, 008(^2/2)}. Set now the initial state of 
subsystem B as a mixed state: Pb — \{\s){s\ + \t^){t+ \ + 
\to) {to I + \t-){t-\), and apply the MPEA to this system. Af- 
ter a series of partial measurements in the basis of |(^^) = 1 1), 
the state of subsystem B evolves to a singlet state |s), which 
corresponds to the largest eigenvalue of Vb{t). We resolve 
the corresponding phase as zero, thus the eigenvalue is one. 
The survival probability, P„i, of the state \ipj^) = |1) on sub- 
system A after m successful measurements, and the fidelity, 
F{m), for subsystem B to be in state |umax) are shown in 
Fig. 2. Since F{m) is close to one, the success probability is 
determined by P,„ . 

If we start from a pure initial state of subsystem B and ap- 
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MPEA, we implement the unitary transformation of the whole 
system only once; the non-Hermitian matrix is constructed as 
the evolution matrix on the target register By performing peri- 
odic partial measurements on the interacting register, the state 
of the target register is driven automatically to an eigenstate 
of the transformation matrix. Using single-qubit tomography 
and mQFT, we can obtain the complex eigenvalues of the non- 
Hermitian matrix. This algorithm can be used to study open 
quantum system and in developing other new quantum algo- 
rithms. 



FIG. 2: Survival probability P{m) (•) and fidelity F(rn) (A) for 
\Va) ~ U) for the Jaynes-Cummings model. 



ply MPEA, the state of B can evolve to other eigenstates of 
Vb- For example, applying the MPEA to the above system 
using as the initial state for subsystem B, by perform- 
ing partial measurements in terms of = |1) on A, the 
state of the system B would remain in the \t^) state. We can 
retrieve the real part of the phase factor of the correspond- 
ing eigenvalue up to an accuracy of 2, 8 and 16 binary dig- 
its, respectively, and obtain the eigenvalues as exp[— 0.5177 — 
^27^(0.25 ± 0.25)], exp[-0.5177 - i27r(0.160 ± 0.008)], and 
expj-0.5177-i27r(0.15918±0.00003)], assuming we obtain 
the imaginary part through QST and projective measurement 
exactly. The true eigenvalue is exp(— 0.5177 — i). 

Constructing the Hamiltonianfor performing the controlled 
unitary operation - On the MPEA circuit, we set the control 
qubit as a single spin and label it as subsystem C. Thus, the 
controlled Hamiltonian of the whole system becomes 



H = 



i(l - af;)H =\{l- <y^c){woh'^b + {al + a^) + 
J 



(4) 



This Hamiltonian contains three-body interactions and can- 
not be implemented directly. One could decompose the three- 
body interaction into two-body interactions 1 16, 17] and then 
implement them. Generally, an arbitrary unitary evolution 



matrix U 



iHt 



can be decomposed into tensor products 



19(1, which cor- 



of unitary matrices of 4 x 4 and 2x2 1118, 
respond to two- and single-qubit operations respectively, and 
can be implemented on a universal quantum computer 

Conclusion.-We have presented a measurement-based 
quantum phase estimation algorithm to obtain the eigenvalues 
and corresponding eigenvectors of non-Hermitian matrices. In 
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